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ABSTRACT. - Age, growth and mortality of the European pilchard, Sardina pilchardus (Walbaum, 1792) were deter¬ 
mined from 503 specimens collected on the northwestern Aegean coast of Turkey from November 2005 to October 2006. 
Fork length ranged from 8.7 to 14.3 cm while weight varied between 7.05 and 32.66 g. The von Bertalanffy growth equa¬ 
tion was fitted on the basis of mean length-at-age data resulting in parameter values of Lo= = 15.23 cm, K = 0.47, t„ = -1.21. 
The b value for males ( b = 2.455) was somewhat lower than that for females ( b = 2.642), showing negative allometry. The 
oldest male and female were estimated to be three years old. Total (Z), natural (M) mortality, fishing (F) mortality were 
0.94,0.64,0.30 per year, respectively. The exploitation ratio E = 0.32 indicated that the population was slightly exploited. 
The sex ratio was skewed in favour of females (1:1.31). The monthly values of gonadosomatic index (GSI) of females indi¬ 
cated that spawning occurred mainly between September and May with a peak in February. 


RESUME. - Age, croissance et mortalite de la sardine, Sardina pilchardus dans la baie d’Edremit (Turquie). 

L’age, la croissance et la mortalite de la sardine, Sardina pilchardus (Walbaum, 1972) ont ete determines a partir 
d’echantillons (N = 503) collectes sur les cotes turques du nord-ouest de la mer Egee, entre novembre 2005 et octobre 
2006. La longueur a la fourche a varie entre 8,7 et 14,3 cm, et le poids entre 7,05 et 32,66 g. L’equation de croissance de 
von Bertalanffy, basee sur les parametres de longueur moyenne par age, a donne les resultats de = 15,23 cm, K = 0,47 
et t Q = -1,21. Le poids s’accroit de fajon allometrique pour la femelle et le male ( b = 2,642 et b = 2,455, respectivement). 
L’age des plus vieux individus, males et femelles, a ete estime a trois ans. Les mortalites totale (Z), naturelle (M) et de 
peche (F) ont ete de 0,94,0,64 et 0,30 par an, respectivement. Le taux d’exploitation E = 0,32 indique que la population 
n’est que faiblement exploitee. La sexe ratio est en faveur des femelles (1 : 1.31). Les valeurs mensuelles de l’indice gona- 
dosomatique (GSI) des femelles suggerent que la ponte a lieu entre septembre et mai, avec une pointe en fevrier. 
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The sardine, Sardina pilchardus (Walbaum, 1792) is 
a pelagic fish species whose distribution extends generally 
from southern Morocco to British Channel, and from the 
Azores to the eastern Mediterranean and Black Sea (Parrish 
et al., 1989). 

The European sardine is one of the most important small¬ 
sized pelagic species in Turkish waters. It has high com¬ 
mercial importance, being targeted by purse-seine fisheries 
across most of its distribution area and by pelagic trawlers 
mainly in Turkish waters. Namely, the mean annual catch 
of sardines raised up from 12,000 to 20,941 tonnes between 
2003 and 2007, as main target of beach seine on the Turkey 
coasts (Turkstat, 2008). Sardine is an important species, 
widely caught: 15,681 tonnes in Edremit Bay during years 
1995-1999 (Kara and Gurbet, 1999). 

In Europe, Quintanilla and Perez (2000) reported that 
sardina population is mostly composed of age group 2, due 
to changes in the distribution pattern by age. Cabral et al. 
(2003) investigated the abundance of sardine as the main 
target species in the central coast of Portugal. Silva (2003) 
gave morphometric variation among sardine populations 


from the northeastern Atlantic and the western Mediterra¬ 
nean with some biological characteristics. Voulgaridou and 
Stergiou (2003) discussed the decline in the mean TL and 
the increase in the slope b of the length-weight relation¬ 
ship, which are in relation to the recent increase in the effort 
expended on sardine populations and within the framework 
of density-dependent effects on growth in the Eastern Medi¬ 
terranean Sea. Somarakis et al. (2006) reported that the peak 
of the spawning period occurred rather in the Aegean Sea 
(December) than in the less-productive Ionian Sea (Febru¬ 
ary). Sinovcic et al. (2008) investigated the population struc¬ 
ture of sardine stock from the Krka River estuary. 

As to Turkish Seas, there are some scarce published stud¬ 
ies concerning various aspects of sardine biology, ecology, 
genetics, and fisheries (reproduction: Cihangir, 1996; Akyol 
etal., 1996; age and growth: Mater and Bayhan, 1999; Kara- 
kayis and Togulga, 2000; and stock structure: Sarmasik et 
al., 2008). Nevertheless, there is an absence of regular and 
accurate data on growth, maturity, mortality of sardine in 
Edremit Bay, northern Aegean Sea. 

With regard to commercial importance of sardine popu- 
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lation in Turkish commercial fishery, particular attention has 
been given to biological aspects for a proper utilization and 
management purposes. The aim of this paper was to examine 
the population structure of sardine such as growth, length- 
weight relationship and condition, in order to provide better 
knowledge and comparisons with the relevant studies, and 
thus improve protection of the sardine stock in the Edremit 
Bay, northern Aegean Sea. 

MATERIALS AND METHODS 

Fish were obtained from commercial purse seining 
between November 2005 and October 2006 in Edremit Bay, 
northern Aegean Sea (26°57’-26°34’E and 39"17'-39°34' N) 
(Fig. 1). A total of 503 specimens were obtained from dif¬ 
ferent boats each month and about 50 samples were selected 
randomly. Specimens were measured to the nearest 0.1 cm 
fork length (FL) and weighed to the nearest 0.1 g. 

Sagittal otoliths were removed and kept in dry enve¬ 
lopes. Age was read from the whole otoliths immersed in 
glycerin (25%) and alcohol (75%). As a rule, only right oto¬ 
liths were used for age determination. The number of opaque 
zones (summer rings, appearing bright under reflected 
light) and the presence of marginal translucent zone (win¬ 
ter rings, appearing dark under reflected light) were checked 
by two readers. In order to avoid subjectivity effect on age 
estimation, there was an interval of 1 month between read¬ 
ings. Translucent bands that were continuous around the 
entire otolith circumference were considered annuli and 
the total number of these bands was recorded as age. Age- 
classes were assigned based on the number of annuli and the 
month the fish were collected (Chugunova, 1963; Gordoa et 
al., 2000). Length-weight relationships were calculated by 
applying an exponential regression equation W = aU\ where 
W is the weight (g), L is the fork length (cm), and a and b are 
constants. Von Bertalanffy growth equations were calculated 
according to L t = Loo[l - e~ k(t -to)] for FL, where L t is the 


length of fish in cm at age t, Loo is the asymptotic fish length 
in cm, e is the base of natural log (2.71828), t is the fish age 
(year), to is the hypothetical time at which the length of the 
fish is zero, K is the rate at which the growth curve reaches 
the asymptote (Ricker, 1975; Sparre and Venema, 1992). 

Fulton’s condition factor (CF) was calculated as 
CF = [W / L 3 ) x 100 for each sex to assess the maturity and 
condition of specimens (Dulcic et al., 2000). 

The index of growth performance was calculat¬ 
ed according to Pauly and Munro (1984) as follows: 
<t>= Logiok + 21ogioLoo. where, k and Loo are Von Bertalanffy 
growth equation parameters. 

Sex was determined by macroscopic observation of the 
gonads. The overall sex ratio was determined. Deviations 
from 1:1 hypothesis were statistically tested by chi-squared 
analysis (Sokal and Rohlf, 1981). Spawning period was 
determined by analysing the monthly percentages of mature 
individuals (on the basis of macroscopic classification) and 
the mean gonadosomatic index (GSI) over the one-year 
study period. Individual GSI was calculated by the equation: 
GSI = (GW / SW) x 100, where GSI is the gonadosomatic 



Sea. 



Figure 2. - Fork length and weight frequency distribution of 5. pilchardus. 
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Table I. - Mean fork length (FL, cm), and mean weight (W, g), mean condition factor (CF), standard error (SE) for different age groups of 
S. pilchardus males and females. Min-max = minimum and maximum values; N = number of fish; P-value: significance level. * Signifi¬ 
cance values belong to length and weight values. ** Significance values belong to CF values. 





Males 




Females 




Age 

N 

FL ± SE 

(min-max) 

W ± SE 

(min-max) 

CF ± SE 

(min-max) 

N 

FL ± SE 

(min-max) 

W ± SE 

(min-max) 

CF ± SE 

(min-max) 

p = 0.05* 

p = 0.05** 

I 

12 

98.67 ±2.12 

7.18 ±1.42 

1.19 ±0.09 

13 

10.0 ± 1.24 

11.51 ±0.78 

1.11 ±0.056 

p >0.05 

p >0.05 

II 

183 

(8.7-10.4) 

11.91 ±0.43 

(7.18-23.64) 

19.34 ±0.21 

(0.89-2.10) 

1.15 ±0.01 

239 

(8.7-10.9) 

12.01 ±0.39 

(7.05-16.12) 

19.89 ±0.22 

(0.79-156) 

1.14 ±0.01 

p >0.05 

p>0.05 

III 

23 

(10.6-13.5) 

13.22 ±0.70 

(10.74-32.66) 

24.65 ±0.61 

(0.67-2.01) 

1.08 ±0.02 

33 

(10.5-13.5) 

13.19 ±0.81 

(10.76-29.9) 

24.76 ±0.63 

(0.65-2.01) 

1.08 ±0.02 

p >0.05 

p>0.05 


218 

(12.7-14.3) 

(17.94-30.83) 

(0.78-1.28) 

285 

(12.7-14.3) 

(21.63-32.66) 

(0.89-1.28) 




Table II. - Age-length key for S. pilchardus in Edremit Bay, based on oto¬ 
lith readings. 


TL (cm) 

Age groups 

Total 

I 

II 

III 

8.7 

3 



1 

9.0 

2 



4 

9.5 

1 



21 

10.0 

14 



66 

10.5 

5 

27 


80 

11.0 


52 


100 

11.5 


131 


90 

12.0 


116 


53 

12.5 


89 

7 

28 

13.0 


5 

44 

9 

13.5 


2 

3 

5 

14.0 



2 

2 

TOTAL 

25 

422 

56 

503 

% N 

4.97 

83.89 

11.13 

100 

FL 

99.8 ± 1.28 

11.9 ±0.28 

13.20 ±0.40 

12.00 ±0.82 

W 

11.78 ±1.12 

19.65 ±0.24 

24.71 ±0.53 

19.81 ±0.18 

Females 

13 

239 

33 

285 

Males 

12 

183 

23 

218 

F:M 

1:1.08 

1:1.31 

1:1.43 

1:1.31 


index, GW is gonad weight and SW is somatic weight (body 
weight minus gonad weight). Sex ratio was analysed month¬ 
ly (on the basis of macroscopic classification). Deviations 
from 1:1 null hypothesis were statistically tested by chis- 
quared analysis. Total mortality rate (Z) was estimated using 
the following equation (Beverton and Holt, 1957): Z = 1/ 
(t -1’), where t is the average age of the samples and t’ is the 
age at the smallest length of the fish. Natural mortality (M) 


was estimated for shoaling fish using Pauly’s empiri¬ 
cal formula (Pauly, 1980): 

M = 0.8*exp(-0.015 - 0.28LnLoo + 0.65FnK + 0.46FnT°C) 
where Loo and K are the parameters derived from Von 
Bertalanffy equation and T°C the mean annual envi¬ 
ronmental temperature at the surface of the study area 
(10°C). Following the estimation of Z and M, the fish¬ 
ing mortality rate (F) was estimated as: F = Z - M, 
and the exploitation rate as E = F / Z (Pauly, 1980). 

RESULTS 

Length and weight frequency distribution 

The fork length of all individuals (n = 503) ranged 
from 8.7 to 14.3 cm and weight from 7.05 and 32.66 g 
(Fig. 2). Mean fork length and weight were similar 
for females and males (Tab. I), and the differences 
were not significant (t-test, p >0.05). 


Age composition and sex ratio 

Age and sex distribution data are summarized in 
table I. Age of captured fish ranged between I and 
III, while the second year classes were dominant. 
Because of selectivity of the nets, the ‘0’ age group 
was not represented in the samples. There were about 
43% males and 57% females, and differences between 
sexes according to age were not statistically signifi¬ 
cant (p >0.05). The sex ratio (F:M) was 1:1.31. 

Growth 

According to von Bertalanffy growth equations of all 
individuals, a theoretically maximal length of 15.23 cm is 
realistic because the largest specimen sampled during the 
surveys was 14.3 cm. The age-length key is given in table II. 
Age classes covered in the study by otolith sample analysis 
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Figure 3. - Length- weight relationships in females (A) and males 
(B) of S. pilchardus. 


ranged from I to III. Growth performance index (<t>) of Sar¬ 
dina pilchardus is estimated as 2.04. 

Length-weight relationships 

The fork length-weight relationships were evaluated 
separately for males and females (Fig. 3). The calculated 
length-weight equation for females was: W = 0.0274L 2 - 642 
(r 2 = 0.6269) and for males, W = 0.0435L 2455 (r 2 = 0.6396). 
Weight increased negative allome trie ally with size since the 
value of b = 2.642, b = 2.455, respectively had a significant 
difference from the value 3.0 (t-test, p < 0.05) (Fig. 3). 

Condition factor 

The mean condition factor for females was similar to that 
of males, and the differences between sexes were not signifi¬ 
cant (p >0.05, t-test. Tab. I). 

Gonadosomatic index 

Gonad development was followed using the GSI%. 
Monthly changes are plotted in figure 4. Spawning occurred 
between September and May with a peak in February. Dif¬ 
ferences between values of GSI according to months, espe¬ 
cially April 2006, are statistically significant (d.f.: 7.836; F: 
43.252; p < 0.05; one-way ANOVA). 

Mortality 

Total mortality for combined sexes was Z = 0.94 year' 1 . 
The natural mortality was found to be M = 0.64 year 1 . Then, 
the calculation of fishing mortality gave F = 0.30 year 1 . The 
exploitation rate was computed as E = 0.32. 

Discussion 

In this study, a total of 503 specimens of S. pilchardus 
from the Edremit Bay were examined from November 2005 



Nov. Dec. Jan. Feb. Mar. Apr. May Sep. Oct. 

2005 2005 2006 2006 2006 2006 2006 2006 2006 

Months 

Figure 4. - Average of monthly gonadosomatic index (Ave. GSI) of 
female S. pilchardus, from Edremit Bay. 

to October 2006. Sardine is a short living species in Edremit 
Bay. Sardine size ranged from 8.7 to 14.3 cm and the age of 
captured fish between groups I and III, having young indi¬ 
viduals. These ranges vary with different studies, apart from 
the similar to those recorded in Aegean Sea, which is close 
to our studying area (Mater and Bayhan, 1999; Karakayis 
and Togulga, 2000; Silva et al., 2008) (Tab. III). This vari¬ 
ation may be due to different stages in ontogenetic develop¬ 
ment, as well as differences in condition, length, age, sex, 
gear selectivity, gonadal development, and geographical var¬ 
iations (Ricker, 1975; Tirasin, 1993). In addition, as seen in 
table III, the relevant literature reported the ages to IV. But, 
the fact that fishing by purse seine has not been forbidden 
in Edremit Bay has lead to catch young sardines named as 
papalina, which are chosen as target species (Ceyhan et al., 
2006). 

As seen, the population consisted of 57% females and 
43% males. Although sex ratio in most species was close to 
1, this may vary from one species to another, differing from 
one population to another of the same species, and may also 
vary year after year within the same population (Nikolsky, 
1963). While sex ratio in sardine was found being close to 
1 in some studies (Mater and Bayhan, 1999), it was pointed 
out that the sex ratio of females was higher than males (Kar¬ 
akayis and Togulga, 2000). 

The exponent b of the length-weight relationships in 
each sex (b = 2.642 for females, b = 2.455 for males) showed 
that weight increased with length in negative allometry. b 
values in the Edremit Bay were found to be close to esti¬ 
mates in Adriatic and Ionian Seas (Petraikis and Stergio, 
1995; Koutrakis and Tsikliras, 2003; Mustac and Sinovcic, 
2007), but also different from the data found in other stud¬ 
ies with positive allometry in Tunisian waters by Khemiri 
(2006) (Tab. III). The b values are often 3.0 and generally 
comprised between 2.5 and 3.5. As fish grows, changes in 
weight are relatively greater than changes in length, due to 
approximately cubic relationships between fish length and 
weight. The b values in fish differ according to species, 
sex, age, seasons and feeding (Ricker, 1975; Bagenal and 
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Table III. - Age structure, parameters of length-weight relationships (a and b), growth (L™, K, to, < E > ) and CF of S.pilchardus in this and previous studies. - indicates absence of 
data. 
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Table IV. - Spawning seasons of S.pilchardus at various localities and average temperatures according to previous studies. 


Months 


References 
Re etal. (1990) 

Hossucu (1992) 

Sola etal. (1992) 

Cihangir, 1996 

Karakayis and Togulga (2000) 
Guisande et al. (2001) 

Zwolonski et al. (2001) 

Tekinay et al. (2002) 

Ettahiri et al. (2003) 

Voulgaridou and Stergiou (2003) 
Stratoudakis et al. (2004) 

Bellier et al. (2004) 

Tsikliras and Antonopoulou (2006) 
Silva et al. (2006) 

Amenzoui et al. (2006) 

This study (2005-2006) 


J 

F 

M 

A 

M 

J 

J 

A 

S 

o 

N 

- 

- 

- 

- 





- 

- 

- 


D 


Locality 
Portugal 
Izmir Bay 
Cantabrian Sea 
Aegean Sea 
Izmir Bay 
Atlantic-Iber. 
Portugal 
Dardanels 
Northern Atlantic 
Thermaikos Gulf 
Northern Atlantic 
Biscay Bay 
NorthemMedit.Sea 
NorthemAtl antic 
Morocco 
Edremit Bay 


Table V. - Total (Z), natural (M) and fishing (F) mortality rate and 
the exploitation ratio (E) of S. pilchardus in this and previous stud¬ 
ies. 


Literature 

Z 

M 

F 

E 

Zupanovic (1955) 



0.17 


Sinovcic (1986) 


0.5 



Sinovcic (1991) 


0.3 



Pertierra and Perrotta (1993) 


0.29-0.62 



Voulgaridou and Stergiou (2003) 

5.65 

1.34 

4.31 

0.77 

Santoyanni et al. (2005) 


0.5 



This study (2006) 

0.94 

0.64 

0.30 

0.32 


Tesch, 1978). In addition, changes in fish shape, physiologi¬ 
cal conditions, different amounts of food available, life span 
or growth increments can all affect the b growth exponent 
(Frost, 1945; Treer et al., 1998, 1999; Koc et al., 2007). 
Changes in fish shape, physiological changes, hydrological 
environmental conditions, different food availability during 
life span, growth increment can all affect growth exponent 
b (Sinovcic, 2000). The b values in fish differ according to 
species, sex, age, seasons and feeding (Ricker, 1975; Bage- 
nal and Tesch, 1978). 

The theoretical maximal length of 15.23 is realistic 
because the largest specimen sampled during the survey was 
14.3 cm. The maximal recorded lengths prior to our study 
are given in table III. Weight-at-age estimates were more 
variable as a measure of growth length estimates. 

Geographic location and environmental conditions such 
as temperature, organic matter, quality of food, time of cap¬ 
ture, stomach fullness, disease, parasitic loads (Bagenal and 
Tesch, 1978), temperature, organic matter, quality of food 


and the water system in which the fish live (Wootton, 1992; 
Treer et al., 1998, 1999) can also effect weight-at-age esti¬ 
mates. 

In our estimations, the asymptotic fork length of sardine 
was L<x = 15.23 and K value was 0.47 for all individuals, 
indicating that sardine is a fast growing species in Edremit 
Bay. The estimations of L y _ and K values in relevant litera¬ 
ture are shown in table III. A trade-off between growth rate 
(K) and maximum size (Loo) is often found. The trade off is 
influenced by several factors, like temperature, mortality or 
food availability. Increased food availability causes a shift 
towards larger maximum size, but may not increase growth 
rate (Tserpes and Tsimenides, 2001). Growth perform¬ 
ance index of Sardina pilchardus estimated in this study 
(ct» = 2.04) confirms the relevant literature (Tab. III). 

Our findings in CF are higher than those estimated in the 
literature (Tab. III). Edremit Bay, in the northern Aegean 
Sea, is a place where two currents meet and is rich in plank¬ 
ton because of upwelling. In addition, the area is fed by 
waters rich in nutrient from erosion through the Sea of Mar¬ 
mara and the Black Sea. With these facilities, Edremit Bay 
has a high potential for sea food, especially fish (Bilecik, 
1989; Togulga, 1997). It can be pointed out that sardine in 
the Edremit Bay was in good condition. 

In the present study, the mean monthly gonadosomatic 
index values were the highest in February for Edremit Bay 
population (Fig. 4). Spawning occurred between September 
and May. Spawning seasons in Edremit Bay are in accord¬ 
ance with relevant studies, except for those of Guidande 
et al. (2001) in Atlantic-Iberian Peninsula. Due to differ¬ 
ent ecological and climate conditions, the starting and fin- 
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ishing time of reproduction may include different months. 
Spawning periods of fish vary with respect to their species; 
the ecological characteristics of fish are determined by such 
ecological differences as stagnant or running water, as well 
as altitude, temperature and quality of food (Nikolsky, 1963) 
(Tab. IV). 

The high exploitation rate (E = 0.76) estimated for the 
study period in the northwest Aegean Sea, indicated that 
intensive fishing, by removing the largest individuals, may 
be implicated in such a decline in maximal and mean TL 
of sardine during 1996-2000 (Voulgaridou and Stergiou, 
2003) (Tab. V). The exploitation rate for the studied period 
(E = 0.30), which is lower than the expected optimal exploi¬ 
tation level (E = 0.50), revealed light to moderate exploita¬ 
tion of stocks in the studied area. 

As a conclusion, our results lead to consider that the effi¬ 
cient management of the fisheries of small pelagic species 
along Edremit Bay would be obtained by creating marine 
protected areas on spawning grounds and in directing fish¬ 
ing efforts farther offshore. The fisheries strategy should be 
planned so that the fishing period follows the reproductive 
period. 
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